Phase diagram of hard colloidal platelets: a theoretical account 
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We construct the complete liquid crystal phase diagram of hard plate-like cylinders for variable 
aspect ratio using Onsager's second virial theory with the Parsons-Lee decoupling approximation 
to account for higher-body interactions in the isotropic and nematic fluid phases. The stability of 
the solid (columnar) state at high packing fraction is included by invoking a simple equation of 
state based on a Lennard- Jones-Devonshire (LJD) cell model which has proven to be quantitatively 
reliable over a large range of packing fractions. By employing an asymptotic analysis based on 
the Gaussian approximation we are able to show that the nematic-columnar transition is universal 
and independent of particle shape. The predicted phase diagram is in qualitative agreement with 
simulation results. 
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I. INTRODUCTION 

Many colloidal dispersions, such as natural clays, and 
macromolecular systems consist of oblate or disk-shaped 
mesogens whose intrinsic ability to form liquid crystalline 
order gives rise to unique rheological and optical prop- 
erties. Despite their abundance in nature, the statistical 
mechanics of fluids containing non-anisometric particles 
in general (and oblate ones in particular) has received far 
less attention than that of their spherical counterparts. 
The possibility of a first order disorder-order transition 
from an isotropic to a discotic nematic fluid of platelets 
was first established theoretically by Onsager [l| in the 
late 1940s. Although originally devised for rod-like parti- 
cles in solution, his theory also makes qualitative predic- 
tions for plate-like particles based on the central idea that 
orientation-dependent steep repulsive interactions alone 
are responsible for stabilising nematic order. 

The intrinsic difficulty with platelets, as pointed out 
by Onsager in his original paper, is that the contribu- 
tion of third and higher body correlations can no longer 
be neglected like for thin rod- like species. Consequently, 
the original second-virial treatment is expected to give 
qualitative results at best Q. In a pioneering simula- 
tion study, Eppenga and Frenkel 0] provided numerical 
evidence for an isotropic-nematic transition in systems 
of infinitely thin circular disks and found the transition 
densities to be much smaller and the first order nature of 
the transition to be much weaker than predicted by the 
Onsager theory. Owing to the simplicity of the model, 
the discrepancy can be attributed entirely to the neglect 
of third and higher body virial terms in the theory. 

At high densities, the nematic phase of disk-shaped 
particles becomes unstable with respect to columnar or- 
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der, characterised by a planar (2D) hexagonal arrange- 
ment of columns each with a liquid internal structure. 
Similar to the formation of the nematic phase, the sta- 
bility of the columnarphase can be explained solely from 
entropic grounds [1, If. Although the system loses con- 
figurational entropy because of the partial crystallisation 
associated with columnar order, this loss is more than off- 
set by a simultaneous increase in translational entropy, 
i.e. the average available space for each particle increases. 

Attempts to improve Onsager's second virial theory 
have met with variable success (see Ref. @ for a recent 
overview). These approaches usually involve integral 
equation or geometric density functional methods whose 
applicability is often restricted to isotropic fluids @, HI , 
models with parallel or restricted orientations Q or par- 
ticles with vanishing thickness 0, |Tl| ■ A recent gener- 
alisation of the fundamental measure approach towards 
arbitrarily shaped hard convex bodies provides a poten- 
tially promising avenue to address more realistic models 
for liquid crystal ordering [l2T | . The influence of higher- 
body correlations can only be assessed numerically via 
computer simulation [l3j]. For cut spheres, the virial co- 
efficients have been quantified up to the 8th order both 
in the fluid isotropic |14| and nematic state [l5| . Despite 
the large number of virial terms, the convergence of the 
virial expansion of the free energy was found to be insuffi- 
cient to provide an accurate description of dense nematic 
and columnar states. Alternatively, Scaled Particle The- 
ory (SPT) can be used to incorporate higher virial terms 
in an indirect manner. Whilst SPT produces reasonable 
results for infinitely thin disks [l6| , its extension to finite 
aspect ratios leads to poor predictions for the isotropic- 
nematic transition densities [l7l ]. 

A simpler strategy to account for higher-body parti- 
cle correlations in the isotropic and nematic fluid state 
is provided by the so-called Parsons-Lee decoupling ap- 
proximation [l8l . [lqI . I20I ]. The basic assumption of this 
approach is that the pair correlation function g(r) of a 
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fluid of hard anisometric bodies, which depends rather 
intractably on the centre-of-mass distance vector Ar and 
oricntational unit vectors u and u', can be mapped onto 
that of a hard sphere fluid with the same packing fraction 
via: 



g(Ar/a ; u, u') = c/Hs(Ar/cr(Af ; u, u')) 



(1) 



with do some reference distance (e.g. particle diameter) 
and <r(Af ; u, u') the distance of closest approach of a 
pair of hard anisometric bodies at a given set of orien- 
tation unit vectors. In case of hard spheres the distance 
of closest approach is simply the hard sphere diameter 
do. Eq. ([T|) provides a natural route of decoupling the 
translational and oricntational degrees of freedom. Start- 
ing from the generalised virial equation it is possible to 
derive an expression for the excess free energy which is 
similar to the one from Onsager with the particle density 
replaced by a reseated density involving the hard sphere 
excess free energy. Whilst the decoupling approximation 
is known to work well for short hard spherocylinders [2l| , 
its merits for plate-like cylinders have not been investi- 
gated so far. This we intend to do in the present paper. 

As for the columnar state, the high degree of posi- 
tional and oricntational order can be exploited to devise 
simple free- volume approaches inspired by the Lennard- 
Jones-Devonshire (LJD) cell model [HlpJIl]. This was 
first done by Taylor and Hentschke [25|, [26| for the high- 
density liquid crystal states of parallel cylinders which do 
not exhibit an isotropic phase. The approach was further 
developed and modified in Ref . [13 showing that a quanti- 
tatively reliable equation of state for the columnar phase 
can be obtained by accounting for the oricntational en- 
tropy of the particles, neglected in the original version. 
In this paper, we will combine the Onsager-Parsons ap- 
proach for the isotropic and nematic fluid state with the 
modified LJD cell theory for the columnar phase to trace 
the complete phase diagram for freely rotating hard cylin- 
ders as a function of thickness-to-diameter ratio. The 
theoretical predictions will be tested against simulation 
results for hard cut-spheres. In view of the inherent dif- 
ficulty of capturing multi-particle correlations in dense 
plate fluids, the overall performance of the present the- 
ory must be deemed satisfactory. Although quantitative 
agreement with simulation data is generally lacking, the 
theory does manage to reproduce the generic features 
of the phase diagram and provides a simple theoreti- 
cal underpinning for the relative stability of nematic and 
columnar order as a function of the plate aspect ratio. 

This paper is constructed as follows. Section II and 
III are devoted to a detailed exposition of the Onsager- 
Parsons and modified LJD theories, respectively. The 
phase diagram emerging from the present theory will be 
presented and discussed in Section IV. Next, algebraic 
forms of the nematic and columnar free energy are given 
which allow us to obtain universal scaling results for the 
nematic-columnar transition. Finally, some concluding 
remarks are formulated in Section VI. 



II. ONSAGER-PARSONS THEORY FOR THE 
ISOTROPIC-NEMATIC TRANSITION 

Let us consider a system of hard cylinders with length 
L and diameter D in a macroscopic volume V . For disk- 
like cylinders we consider here the aspect ratio L/D is 
much smaller than unity. The particle concentration is 
expressed in dimensionless form via c = ND 3 /V. Fol- 
lowing Ref. H^, the Helmholtz free energy within the 
Onsager-Parsons-Lee approach takes the following form: 



N 



l + (ln47r/(u)) 



2 D 3 



(2) 



with = fc^T the thermal energy (fc^ repre- 

sents Boltzmann's constant and T temperature) and 
V = V / D 3 the dimensionless thermal volume of a 
platelet including contributions from the rotational mo- 
menta. The brackets ((•)) = / du/(u)(-), (((■))) = 
JJ dudu'/(u)/(u')(-) denote single and double orien- 
tational averages involving some unknown distribution 
/(u) of the orientation unit vector u of the plate normal 
which is normalised according to J du/(u) = 1. Several 
cntropic contributions can be distinguished in Eq. 
The first two are exact and denote the ideal transla- 
tional and oricntational entropy, respectively. The last 
term represents the excess translational or packing en- 
tropy which accounts for the particle-particle interactions 
on the approximate level of pair-interactions. The key 
quantity here is the excluded volume V exc i between two 
plate-like cylinders at fixed inter-particle angle 7 with 
sin 7 = |u x u'|. This quantity has been calculated in 



qua 

closed form by Onsager [l| and reads: 



D 3 



TT L /7T 7T 

-| sin 7 1 + — - + 2£7(sin7) + - cos 7 



-2|* 
D 



D \2 



sin 7 1 



(3) 



with E{x) the complete elliptic integral of the second 
kind. Although the structure of Eq. @ is similar to 
the classic Onsager second- virial free energy, the effect of 
higher order virial terms are incorporated via the scaling 
factor 



G P = 



(1 



(4) 



which depends on the total plate volume fraction = 
c{it/A)L/D. The rescaled density stems from the 
Parsons-Lee method [lH EM H3, HH which involves a 
mapping of the plate pair distribution function onto that 
of a hard sphere system via the virial equation. The 
free energy can ultimately be linked to the Carnahan- 
Starling expression for hard spheres, which provides a 
simple strategy to account for the effect of higher-body 
particle interactions, albeit in an implicit and approxi- 
mate manner. As required, Gp approaches unity in the 
limit 4> —> in which case the original sccond-virial the- 
ory is recovered. 
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Let us now specify the oricntational averaging. By 
definition, all orientations are equally probable in the 
isotropic (I) phase and fj = l/4ir. The oricntational 
entropy then vanishes: 



(ln47r/j)j = 



(5) 



/3tt 


7T^ 




TT 


(i) 


T" 




D 


h 2 





If we use the random isotropic averages ((sin 7))/ = 7r/4, 
((E(sm^f)))i = 7r 2 /8 and ((cos 7})/ = 1/2 the excluded 
volume entropy reduces to: 

((W7)))/ _ 

D 3 8 

With this, the free energy of the isotropic phase is fully 
specified. 

In the nematic (N) phase, the particles on average 
point along a common nematic director n and the orien- 
tation distribution /(u-n) is no longer a trivial constant. 
For a uniaxial nematic phase, /(u) = f(9) involving the 
polar angle < 6 < tt between the plate normal and the 
director, with / being a peaked function around 9 = 
and 9 = it. 

The equilibrium form follows from the minimum con- 
dition of the free energy: 



(7) 



where the Lagrange multiplier A ensures the normalisa- 
tion of /. Applying the condition to Eq. |(5J) leads to a 
self-consistency equation for /: 



m 



exp [~cG P ((j)) j du'D 



^ cxcl (|uxu'|)/(u')] 



Jduexp [-cG P (<f>) J du'D- 3 V excl (\u x fi'|)/(u')]] 

(8) 

which needs to be solved numerically for a given particle 
concentration (29|. Note that the isotropic distribution 
/ = 1/4-7T is a trivial solution of Eq. ([5]), irrespective of 
c. At higher densities, nematic solutions of Eq. (|SJ) will 
appear which give rise to a lower free energy than the 
isotropic one. The nematic order parameter S 1 , defined 
as: 



S = / duV 2 (cos9)f(u) 



(9) 



[where V-zix) = (3x 2 — l)/2] is used to distinguish the 
isotropic state (S = 0) from the nematic (0 < S < 1). 
Once the equilibrium oricntational distribution function 
is known, the pressure P = —(dF/dV)N.T and chemical 
potential /i = (dF/dN)v,T can be specified to establish 
phase equilibria between isotropic and nematic states. 



III. LJD CELL THEORY FOR THE 
COLUMNAR PHASE 

To describe the thermodynamical properties of a 
columnar phase we use an extended cell theory as pro- 
posed in Ref. [27]. In this approach, the structure of a 



columnar phase is envisioned in terms of columns ordered 
along a perfect lattice in two lateral dimensions with 
a strictly one-dimensional fluid behaviour of the con- 
stituents in the remaining direction along the columns. 
As for the latter, the configurational integral of a system 
of parallel platelets with thicknesses L and diameter D 
with their centre-of-mass moving along the plate normal 
on a line of length I is formally given by [30[ : 



Q fluid (A,£,T) 



1 



K N N\ 



[i - NL] 



N 



(10) 



with A the thermal de Broglie wavelength. The columns 
are assumed to be strictly linear and rigid so that fluc- 
tuations associated with bending of the columns can be 
ignored. Next, we allow the platelets to rotate slightly 
about their centre-of-mass. At high packing fractions, 
the rotational freedom of each platelet is assumed to 
be asymptotically small and the configurational integral 
above may be approximated as follows: 



hwd(N,e,T) 



Qo 



[£-n(l cS ) 



iN 



(11) 



where Vi represents the total ID thermal volume includ- 
ing contributions arising from the 3D rotational momenta 
of the platelet. Furthermore, Q OI = cxp[— AT(ln47r/)] 
is an oricntational partition integral depending on the 
oricntational probability distribution /. In the mean- 
field description implied by Eq. (fTTj) there is no cou- 
pling between the orientational degrees of freedom of the 
platelets. The rotational freedom of the platelets is ex- 
pressed in an effective entropic thickness, defined as 



(Lea) 



Lil 



1 D 



d(cos6)\6\f(6) 



(12) 



up to leading order in the polar angle 9 which measures 
the deviation of the plate normal from the direction of 
the column. A prefactor of '1/2' in Eq. (fT2|) has been 
included to correct in part for the azimuthal rotational 
freedom and captures the effect that the excluded length 
between two platelets at fixed polar angles becomes min- 
imal when the azimuthal orientations are the same. The 
free energy of the ID fluid then follows from (3F = — In Q: 



[3F a 



aid 



A 



= lnVi/5-l + (ln47r/)-ln l-p{L eS ) (13) 



in terms of Vi = Vi/L, the reduced linear density p = 
NL/i and effective thickness L c g = L e g/L. 

Similar to the nematic case in the previous Section, 
the equilibrium form f(9) is found by a formal minimi- 
sation of the free energy under the normalisation con- 
straint. The corresponding stationarity condition is given 
by Eq. ([7]) . Since the free energy Eq. flB"]) depends only on 
one-particle orientational averages, the equilibrium ODF 
can be obtained in closed form and turns out be of a 
simple exponential form: 



/(0) = |UxpK|0|] 



(14) 
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with 



3D S 

2i7y : 



(15) 



The oricntational averages are now easily carried out and 
the leading order expressions for the orientational en- 
tropy and entropic thickness are given by: 



(In 47T/) = 21n£-2 



(16) 



As a measure for the orientational order along the col- 
umn, we can define a ncmatic order parameter S [of. 
Eq. M] related to £ via: 



S= (T 2 (cos9)) 



(17) 



Let us now turn to the free energy associated with the 
positional order along the lateral directions of the colum- 
nar liquid crystal. A formal way to proceed is to map the 
system onto an ensemble of N disks ordered into a 2D lat- 
tice. Near the close packing density, the configurational 
integral of the system is provided in good approxima- 
tion by the LJD cell theory [H, IH, 0, El- Within the 
framework of the cell model, particles are considered to 
be localised in 'cells' centred on the sites of a fully occu- 
pied lattice (of some prescribed symmetry). Each parti- 
cle experiences a potential energy w"" u (r) generated by its 
nearest neighbours. In the simplest version, the theory 
presupposes each cell to contain one particle moving inde- 
pendently from its neighbours. The A^-particle canonical 
partition function can then be factorised as follows: 



fLJD 



(N) 



-± J dr N eM-PU(r N ) 



^2 / d2rex P 



-~<eii(r) 



N 



(18) 



For hard interactions, the second phase space integral is 
simply the cell free area available to each particle. If we 
assume the nearest neighbours to form a perfect hexago- 
nal cage, the free area is given by Af Tee = \/3(Ac — D) 2 /2 
with Ac the nearest neighbour distance. The configura- 
tional integral then becomes 



(AO 



Afr, 



A 2 



iV 



^A 2 C 
A 2 



JV 



1 



Ac 1 



2N 



(19) 

where the (lateral) spacing Ac = Ac/D is a measure 
for the translational freedom each particle experiences 
within the cage. The free energy associated with the 
LJD cell theory is given by: 



l3F hm 
N 



In 



D 2 



In 



V3 



2 In 



(20) 



The LJD equation of state associated with Eq. ([20]) pro- 
vides a very accurate description of a 2D solid at densities 
near close-packing [32]. If we now apply the condition 
of single-occupancy (i.e. one array of platelets per col- 
umn) we can use Ac to relate the plate volume fraction 
4> = Nvq/V (with vq = (tt/A)LD 2 the particle volume) 
to the reduced linear density p via: 



b*A 2 c = p 



(21) 



with 



in terms of the reduced packing fraction cf>* — Y/ ^ cp 
4>cp = 7r/2"\/3 r* 0.907 the value at close packing. The 
total free energy of the columnar state is now obtained 
by adding the fluid and LJD contributions: 



(3F co 
N 



In Vc - 1 
In'' 1 



2 In 



'A 2 



3D 

Yl \ i 



*A2 



► A 2 



- 21n(l - A^ 1 ) 



(22) 



where the ideal contribution is identical to that of Eq. j2]). 
The final step is to minimise the total free energy with 
respect to Ac- The stationarity condition dF/dAc = 
yields a third-order polynomial whose physical solution 
reads: 



-S^W + 2 1 / 3 A" 2 / 3 

62/3^*^1/3 



with 



K = 27((f>*) 2 + [3(</>*) 3 (32 + 2430*)] 



1/2 



(23) 



(24) 



With this, the free energy for the columnar state is fully 
specified. Unlike the nematic free energy, the columnar 
free energy is entirely algebraic and does not involve any 
implicit minimisation condition to be solved (cf. Eq. 10). 
The pressure and chemical potential can be found in 
the usual way by taking the appropriate derivative of 
Eq. (|2"2"]) . In Sec. [V]we will show that the nematic free 
energy can also be recast in closed algebraic form using a 
simple variational form for the ODF, similar to Eq. (CO 



IV. PHASE DIAGRAM 

Fig. [T] presents an overview of the phase behaviour of 
a hard cylindrical platelets based on the theoretical ap- 
proach described above, along with various simulation 
data for hard cut spheres available in literature. From 
Fig. QJi it is evident that the packing fractions associ- 
ated with the isotropic-nematic coexistence increase for 
larger aspect ratio whereas the nematic-columnar tran- 
sition remains virtually unaffected by the shape of the 
platelet. This observation is in line with the tentative 
phase diagram constructed by Veerman and Frenkel 0]. 
The trends can be understood qualitatively by noting 
that the onset of nematic order occurs if the fraction of 
excluded volume ~ ND 3 /V exceeds a certain universal 
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FIG. 1: (a) Phase diagram for monodisperse colloidal platelets of variable aspect ratio L/D in terms of the plate packing 
fraction (j> = (tt/4)LD 2 N/V. Thin continuous lines serve to guide the simulation data, (b) Dimensionless concentration of 
the coexisting isotropic and nematic phases as a function of the plate aspect ratio. Inset: nematic order parameter S of the 
coexisting nematic phase plotted versus aspect ratio. 



value of about 4 (as reflected in Fig. Q})) whereas colum- 
nar order only becomes stable beyond a critical packing 
fraction, typically cf> ~ 0.4. Whence: 

IN ~ ttL/D, L/D < 1 

0nc - 0.4 (25) 

which implies the presence of a triple aspect ratio, fixed 
by the intersection of both nematic binodals. Although 
at this particular value an isotropic-nematic-columnar 
triphasic coexistence occurs above a certain packing frac- 
tion, the system volume occupied by the nematic phase 
is always infinitesimally small and the situation thus dif- 
fers from a regular tri-phasic coexistence occurring in e.g. 
binary mixtures at a given thermodynamic state point. 
Equating both expression we estimate the triple aspect 
ratio to be L/D = OAir w 0.126, which is very close 
to the value 0.125 obtained from extrapolating the sim- 
ulation binodals from Refs. [33|, S- Beyond the triple 
aspect ratio, the platelets are no longer sufficiently ani- 
somctric to guarantee a stable nematic phase and direct 
transitions from the isotropic fluid to the columnar solid 
occur. We should note that our theory does not take 
into account the theoretically disputed cubatic phase, as 
an intermediate state between the isotropic and colum- 
nar phases. The issue of the stability of cubatic order 
with respect to columnar order is discussed in a recent 
simulation study by Duncan et al. I15J. The transition 
densities from Veerman Q and Bates |35| are systemat- 
ically larger than those reported by Zhang [33[ and van 
der Beek (34J and therefore give rise to a slightly higher 
estimate of the triple value (L/D ~ 0.14). 

The theoretical value L/D = 0.175 deviates consider- 
ably from the ones predicted from simulations, mainly 
because the predicted packing fractions of the coexisting 
nematic and columnar phases are too large. The equa- 
tions of state presented in Fig. [3] demonstrate that the 



main source of error must be the chemical potential of 
cither the nematic or columnar phase nematic branch, 
rather than the pressure. For L/D = 0.05, the pre- 
dicted pressure in the nematic and columnar states are 
fairly close to the simulation results with discrepancies 
less than a few percent in both branches. 

For larger aspect ratios (L/D > 0.1) the agreement be- 
tween theory and simulation is quite satisfactory, despite 
an increased shape difference between the cylinder and 
the cut sphere. For L/D = 0.2, the occurrence of cubatic 
order has been reported in simulation Q which is not 
taken into account in the present model. The isotropic 
binodal point in Fig. QJ, at this value is taken to be the 
mean value between the onset of cubatic order and the 
transition to the columnar state, with the error bar indi- 
cating the boundaries of cubatic order. For the thickest 
species L/D = 0.3, the coexisting high density phase was 
found to be a solid rather than a columnar. Similarly, for 
smaller aspect ratios a continuous columnar-solid transi- 
tion line could be located beyond the nematic-columnar 
moving toward higher packing fraction upon decreasing 
L/D. In our simple cell-fluid model there is, however, 
no distinction between the columnar and solid states due 
to the absence of a freezing transition within the strictly 
ID line fluid representing the structure along the column 
direction. 

The predictive power of the Onsagcr-Parsons theory 
for platelets is perhaps better highlighted in Fig. [TJd, 
where the isotropic-nematic binodals are plotted in terms 
of the reduced number concentration c = ND 3 /V. The 
agreement is reasonable for large aspect ratio but rather 
poor for thin platelets. In the limit of infinitely thin disks 
(L/D — > 0) the coexistence concentrations are identical 
to those obtained from Onsager's second- virial theory viz. 
a = 3.29(16/tt 2 ) and c N = 4.191(16/tt 2 ) [H This is 
easily understood from the fact that the packing fraction 
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FIG. 2: Equation of state for colloidal platelets for two different inverse plate aspect ratios D/L, plotted in terms of the reduced 
pressure PD 3 /k B T versus dimensionless concentration ND 3 /V. (a) Isotropic-nematic density region, (b) Nematic-columnar 
region. 



of infinitely thin disks at a given finite number concentra- 
tion is zero. Consequently, Gp{4>) reduces to unity and 
the Parsons decoupling approximation involving the hard 
sphere excess free energy becomes ineffective. A similar 
reduction to the B% level takes place for infinitely thin 
rods (L/D — > oo). However, contrary to rods, the effect 
of the third virial coefficient B% is finite for disks with 
vanishing thickness. In general, for isotropic systems of 
hard cylinders we have: 



B 3 /B 2 2 
BzlB\ 



0, (L/D -> oo) 

0.444, (D/L -> oo) 



(26) 



where the latter value is taken from Ref. 0. Likewise, 
higher order virial contributions will also be non-zero. 
Simulation studies of the virial terms up to Bj 0, H] 
reveal that higher order virial terms involve alternating 
positive and negative contributions of comparable magni- 
tude, indicating just how complicated virial expansions 
are for dense fluids of platelets. The virial terms gen- 
erated by the Onsager-Parsons free energy can be ob- 
tained from the virial expression for the excess free en- 
ergy (3F CX /N = £„> 2 B„p n - 1 /(n - 1). Applying this to 
Eq. © gives: 



B,, 



(rc + 2)(n-l) /7rX ((W7)» Y 
4(n-2)! \8D D 3 J 



(n > 2) 
(27) 



It is clear that all contributions beyond Bi are zero for 
L/D = thus leading back to the original Onsager result. 
For L/D = 0.1 the reduced third, fourth and fifth virial 
coefficients in the isotropic phase are 0.170, 0.010 and 
3.67 -10 -4 . Comparing these with the numeri cally exact 
values 0.508, 0.111 and -0.10 for cut spheres [lj] shows 
that higher-order correlations are systematically under- 
weighted by the Parsons method. 



ASYMPTOTIC RESULTS FOR THE N-C 
TRANSITION 



A simple rationale for the apparent independence of 
the NC transition with respect to particle shape can be 
obtained by comparing the free energy of the two states 
and exploiting the fact that the nematic order at densities 
close to the transition is very strong. In that case the 
average excluded volume Eq. ([3]) between the particles in 
the nematic phase can be approximated by retaining the 
leading order contribution for small inter-rod angles 7: 



L 



D 



22/ 



«Kxcl(7)» ~ 27T- + - + -5- «?» ( 28 ) 



D 2 



The orientational averages indicated by the brackets can 
be estimated using a Gaussian Ansatz for the ODF [37j : 



((•)> 



d0|sin#| dyf G (e)(-) (29) 
it/2 Jo 



in terms of the following one-parameter Gaussian varia- 
tional function: 



f G {0) =Afexp 



-2<*<a) (30) 



where the normalisation factor Af follows from (1) = 1. 
The variational parameter a is required to be much larger 
than unity so that fo is sharply peaked around 8 = 0. 
In that case, the integration over the polar angle 6 in 
Eq. (|29p can be safely extended to ±00 and sin 8 w 6. In 
the asymptotic limit, the normalization constant is given 
by J\f = a/47r and the double orientational average over 
the angle 7 in the nematic phase is found to be [371 ] : 



((7)) 



7T\l/2 



(31) 
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to leading order in a. Similarly, the orientational entropy 
can be approximated by: 



(In47r/ G 



In a — 1 



(32) 



The nematic order parameter [cf. Eq. (|17p] follows from 
S ~ 1 — 3/a. These algebraic results allow the minimi- 
sation of the free energy with respect to a to be carried 
out analytically and leads to a closed expression for the 
nematic free energy. 

Using the asymptotic expressions above, and introduc- 
ing the volume fraction as a density variable the following 
algebraic form for the Onsager-Parsons nematic free en- 
ergy can be produced: 



> (ln0-l) + ( 2 lnf^^0Gp(0)| 



iV 



L 2 



-2 + 4<j>G P ((t>) 



(33) 



which is a sum of the ideal, orientational and excess parts, 
respectively. Similarly, one may derive for the columnar 
free energy: 



^ - (m^-l ) + (2h/ 3 ^ ' 



N 



2 LI- p 



ln^— ^-21n(l-A^) 



(34) 



with p = 4>*A C and combines the ideal, orientational, 
ID fluid and cell contributions, respectively. The only 
explicit shape dependency is the contribution 2 In D/L in 
the orientational part which is identical in both expres- 
sions and therefore does not affect the NC coexistence 
properties. Solving the coexistence conditions gives the 
universal coexistence values 4>m = 0.500 and <j>c = 0.621 
and pressure (3PD 3 = 11.37(D/L). Furthermore, the 
normalised lateral columnar spacing is Ac = 1.075 and 



the equilibrium variational parameters a = 1.232(Z?/L) 2 
pertaining to the nematic order in the nematic phase 
and columnar phases are given by a = 1.276(1?/ L) 2 and 
£ = 5.830(D/L), respectively. 



VI. CONCLUSIONS 

We have combined the Onsager-Parsons theory with a 
simple LJD cell model to address the phase behaviour 
of hard cylindrical platelets with variable aspect ra- 
tio. The theoretical framework provides a simple, yet 
qualitative underpinning for the competitive stability of 
the isotropic, nematic and columnar states, observed in 
Monte Carlo computer simulations. Upon increasing the 
aspect ratio, the window of stability for the nematic 
phase decreases systematically up to a critical value, 
identified as a triple equilibrium. Beyond this value the 
anisomctry of the plates is too small to warrant a stable 
nematic phase and direct transitions from an isotropic 
fluid to a columnar solid occur. It would be intriguing to 
verify whether the stability of cubatic order, suggested 
by the simulations, can be captured with a similar free 
volume concept. Rather than forming a closed packed as- 
sembly of columns, the cubatic phase must be envisioned 
in terms of interacting finite-sized stacks with random 
orientations. This will be explored in a future study. 
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